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Abstract: Wc address the problem of curvature estimation from sampled com- 
pact sets. The main contribution is a stability result: we show that the gaussian, 
mean or anisotropic curvature measures of the offset of a compact set K with 
positive ^-reach can be estimated by the same curvature measures of the offset 
of a compact set K' close to K in the Hausdorff sense. Wc show how these 
curvature measures can be computed for finite unions of balls. The curvature 
measures of the offset of a compact set with positive yu-reach can thus be ap- 
proximated by the curvature measures of the offset of a point-cloud sample. 
These results can also be interpreted as a framework for an effective and robust 
notion of curvature. 
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Stabilite des Mesures de Courbures 



Resume : Ce travail etudie restimation de la courbure d'ensembles compacts 
echantillonnes. Nous montrons que les mesures de courbures moyenne, de Gauss 
et anisotropes des offsets d'un compact K sufRsamment "regulier" , i.e. a /x-reach 
strictement positif, sont stables : elles peuvent etre estimees par les mesures de 
courbures des offsets de tout compact K' suffisamment proche de K pour la 
distance de HausdorfF. Nous explicitons le calcul de ces mesures de courbures 
dans le cas ou K' est un nuage de points dans R^. Ces resultats fournissent en 
particulier un cadre permettant de definir formellement des notions de courbures 
robustes pouvant etre efFectivemcnt calculees. 

Mots-cles : mesures de courbure, cycle normal, inference geometrique, ^- 
reacli, fonction distance, nuage de points 
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1 Introduction 

Motivation We present in this worii a stable notion of curvature. A common 
definition of curvature considers quantities defined pointwisely on a twice dif- 
ferentiable manifold. However, the objects we have to deal with in practice are 
not twice diffcrcntiable: consider the situation where a physical object is known 
through a sufficiently dense point cloud measured on the object boundary. Intu- 
itively, it seems feasible to infer some meaningful information on the curvature 
of the physical object itself. Let us assume that we know both the accuracy of 
the measure, which is an upper bound on the distances between the measured 
points and their closest points on the physical object boundary, and the sam- 
pling density, which is an upper bound on the distance between the points on 
the object boundary and their closest measured sample points. The fact that 
the measure accuracy and sampling density are below a known small value e can 
be expressed by saying that the Hausdorff distance (see the definition in Section 
[2]) between the assumed physical object and the measured point cloud is less 
than e. Even with the guarantee of a small Hausdorff distance, the knowledge 
of the measured point cloud allows many possible shapes for the physical object 
boundary and it is hopeless, without additional assumptions, to infer the usual 
pointwise curvature quantities on the physical object. 

Indeed, a first difficulty is that the geometrical and topological properties 
of a physical object have to be considered at some scale: for example, if one is 
interested in the shape of a ship hull, it may make sense to see it as a smooth 
surface at a large scale (10~^ meters). However, at a finer scale (10"'' meters) 
this same object may appear with many sharp features near the rivets and small 
gaps between assembled sheets of metals. At a still finer scale, one could even 
consider each atom as a separate connected component of the same physical 
object. 

A second difficulty is related to the "pointwise" character of usual curvature 
definition. If the physical object is known up to e in Hausdorff distance, it seems 
impossible to distinguish between the pointwise curvature at two points whose 
distance is of the order of e. Without very strong assumptions on the regularity 
of the unknown object, it is again meaningless to evaluate a pointwise curvature 
quantity. 

In this work, we overcome the first difficulty by considering a scale dependent 
notion of curvature which consists merely, for a compact set K, in looking at 
the curvature of the r-offsets of K {i.e. the set of points at distance less than 
or equal to r from K). This offset, which can be seen as a kind of convolution, 
in the same spirit of similar operators in mathematical morphology or image 
processing, filters out high frequencies features of the object 's boundary. 

We overcome the second difficulty by considering curvature measures instead 
of pointwise curvature. In the case of smooths manifolds, curvature measures 
associate to a subset of the manifold the integral of pointwise curvatures over 
the subset. However the curvature measures are still defined on non smooth 
objects such as convex sets or more generally sets with positive reach 

In practice, in order to state our stability theorem, one has to make some 
assumption on the unknown physical object, namely the positive fi-reach prop- 
erty defined below. To be more precise, our stability theorem still applies to 
objects whose offsets have positive ^-reaches. 
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Some other approaches assume stronger properties such as smoothness or 
positive reach. But is it legitimate to make such assumptions on an unknown 
object? In practice, the only available informations about the physical object 
appears through the physical measures. A distinctive character of our assump- 
tion on the /x-reach of offsets of physical object is that, thanks to the so-called 
critical values separation theorem [B], it can be reliably checked from the mea- 
sured point sample. In this situation our curvature estimations reflect reliably 
the intrinsic properties of the physical object. 

Related previous works. Due to its applications in geometry processing, 
many methods have been suggested that, given a triangulated surface, are able 
to estimate the curvature of an assumed underlying smooth surface (see |32] for 
a survey). Several authors (for example [2]) compute the curvature of a smooth 
polynomial surface approximating locally the triangulated surface. Closely re- 
lated to our work is [STJ [Ml [151 [16] , which study a general definition of curvature 
measure that applies to both smooth surfaces and their approximation by tri- 
angulated surfaces, based on the so-called normal cycle (defined below). The 
proximity of curvature measure is proved using the powerful notion of flat norm 
between the corresponding normal cycles. These ideas are thoroughly reused in 
the present work. 

In [6], in order to address the question of topology determination through 
Hausdorff approximation, the authors have introduced the class of sets with 
positive /i-reach, which can be regarded as a mild regularity condition. In 
particular, this condition does not require smoothness. If a set has a positive 
/i-reach, or at least if it has offsets with positive /i-reach, it is possible to retrieve 
the topology of its offset from the topology of some offsets of a Hausdorff close 
point sample. A stable notion of normal cone has also be defined on this class 
of sets |5]. More recently, the authors have proved [7] that the complement of 
offsets of sets with positive /i-reach have positive reach. Since the normal cycle 
and the associated curvature measures are defined for sets with positive reach 
[22l [24] we may consider normal cycles of offsets of sets with positive /t-reach. 
This paper develops this idea and gives a stability result whose proof uses the 
fact that the boundary of the double offsets of sets with positive /i-reach are 
smooth surfaces. 

Contributions. Our main result states that if the Hausdorff distance between 
two compact sets with positive /t-reach is less than e, then the curvature mea- 
sures of their offsets differs by less that 0{^/e), using an appropriate notion of 
distance between measures. This is then extended through the critical values 
separation theorem to the case where only one set has positive /t-reach, which 
allows to evaluate the curvature measure of an object from a noisy point cloud 
sample. These results improve on the stability results in [311 [131 [HI [H] , which 
were only limited to the approximation of smooth hypersurfaces by homeomor- 
phic triangulated manifolds. In order to provide a concrete algorithm we give 
the formulas that express these curvatures measures on a union of balls (i.e. 
on an offset of the point cloud) . Closest to our work is [S] , which also gives a 
stability result for curvature measures. The main differences are that our result 
also applies to anisotropic curvature measures, whereas [5] is only limited to the 
usual curvature measures. On the other hand, the stability result for curvature 
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measures in derives from a stability result for so-called boundary measures, 
which holds without any assumptions on the underlying compact set, whereas 
ours requires to assume a lower bound on the ^-reach. While the two results 
seem related at first sight, the proof techniques are drastically different. 

Outline. The paper first recalls definitions and properties related to the dis- 
tance function to a compact set and its gradient, the critical function and the 
/x-reach. Then, in Section [31 one recalls classical notions concerning currents 
and one introduces the notions of normal cycle and of curvature measures. In 
Section m we state the stability theorem for sets with positive /x-reaches. We 
then extend this theorem to the case where only one set has positive /i-reach. 
In Section [51 we give the expressions of curvature measures for unions of balls. 
The last section gives the main steps of the proof of the main Theorem. 

2 Definitions and background on distance func- 
tions 

We are using the following notations in the sequel of this paper. Given X C M'*, 
one denotes by X'^ the complement of X, by X its closure and by dX the 
boundary of X. Given A C M'', ch{A) denotes the convex hull of A. 

The distance function dx of a compact subset K of associates to each point 
a; e R'' its distance to K: 

X I— *■ djilx) ~ min (i(.T, y), 

ijeK 

where d{x, y) denotes the euclidean distance between x and y. Conversely, 
this function characterizes completely the compact set K since K = {x € 
W^\dK{x) = 0}. Note that dx is 1-Lipschitz. For a positive number r, we 
denote by Kr the r-offset of K, defined by AV ~ {x \ dxix) < r}. The Haus- 
dorff distance dH{K,K') between two compact sets K and K' in R'^ is the 
minimum number r such that K C K'^. and K' C Kr- It is not difficult to 
check that the Hausdorff distance between two compact sets is the maximum 
difference between the distance functions associated with the compact sets: 

dH{K,K') = sup \dKix)-dK'ix)\ 

Given a compact subset K oiM.'^, the medial axis Ai [K] of K is the set of points 
in \ K that have at least two closest points on K . The infimum distance 
between K and M.{K) is called, according to Federer, the reach of K and is 
denoted reach(K). reach(K) = if -ftT has concave sharp edges or corners. The 
projection map pk that associates to a point x its closest point pk{x) on K is 
thus defined on \M{K). 

A C^'^ function is a function whose first derivative is Lipschitz. A C^'^ 
hypersurface 5 is a (c? — l)-manifold embedded in such that each point 
of S has a neighborhood which is the regular image (that is the image by a 
function whose derivative has maximal rank) by an injective C^'^ function of 
a neighborhood of in M"*^^. Informally, one can say that a C^'^ surface is a 
surface with bounded curvature, which is strictly stronger than and strictly 
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weaker than C^. An embedded compact manifold is C^'^ if and only it has 
positive reach (Federer). 

2.1 The gradient 

The distance function (Ik is not differentiable on /A{K). However, it is possible 
[28] to define a generalized gradient function : K'' — > R"^ that coincides with 
the usual gradient of dx at points where dx is differentiable. For any point 
X GMf^X K, we denote by F^f (x) the set of points in K closest to x (Figure [T|): 

TKix) ^ {y e K \ d{x,y) ^ dK{x)} 

Note that Tk{x) is a non empty compact set. There is a unique smallest closed 
ball aK{x) enclosing Tk{x) (cf. Figure [1]). We denote by 9k{x) the center of 
(T/<-(a;) and by !Fk{x) its radius. 0k{x) can equivalently be defined as the point 
on the convex hull of r/^(.T) nearest to x. For x <E \ K, the generalized 




Figure 1: A 2-dimensional example with 2 closest points, 
gradient V_R-(a;) is defined as follows: 

-^.-eiiix) 



^Kix) 



dxix) 



It is natural to set V/f (x) = for a; e K. Note that for a; G K'' \ K, \\\/k{x)\\ 
is the cosine of the (half) angle of the smallest cone with apex x that contains 
Tk{x). 

2.2 Critical points, critical function and the ^-reach 

The critical points of dx are defined as the points x for which V/f (x) = 0. 
Equivalently, a point x is a critical point if and only if it lies in the convex hull 
of Tk{x). When K is finite, this last definition means that critical points are 
precisely the intersections of Delaunay fc-dimensional simplices with their dual 
{d — /s)-dimensional Voronoi facets [H]. Note that this notion of critical point 



INRIA 



Stability of Curvature Measures 



7 



is the same as the one considered in the setting of non smooth analysis [13| and 
Riemannian geometry [12 . 

The results of this paper rely strongly on the notions of critical function and 
fi-reach, introduced in [5]. 

Definition 1 (critical function) Given a compact set K C M"^, its critical 

function xk ■ (0, +oo) R+ is the real function defined by: 

XK{d)= inf IIVkII 

dK id) 

The function xk is lower semicontinuous. The /x-reach of a compact set K is 
the maximal offset value d for which XK{d') > ^ for d' < d. 

Definition 2 (/i-reach) The /i- reach r^{K) of a compact set K izW^ is de- 
fined by: 

r^iK) = inf{d I XK{d) < fi} 

We have that ri{K) coincides with the reach introduced by Federer [52] • The 
critical function is in sonic sense "stable" with respect to small (measured by 
Hausdorff distance) perturbations of a compact set [Sj . That implies the follow- 
ing theorem [S]: 

Theorem 1 (critical values separation theorem) Let K and K' be two com- 
pact subsets ofW^, dniK, K') < e and ^ be a non-negative number. The distance 
function dx has no critical values in the interval '\Ae/ii^,rf^{K') — 3e[. Besides, 
for any fi' < fi, xk is larger than fi' on the interval 



4e 



— ,r^{K')-3^er^{K') 



2.3 Complement of offsets and double-offset 

It has been proved in that the complement of the offset K,. has positive 
reach for any value < r < r^. Moreover, one has a lower bound for the critical 
function of K^: 

Theorem 2 For r G (0,rp), one has 

reach(K5:) > /ir (1) 

Moreover for any t G (/ir, r), 

2^r~til + ^^) 

Let K C M'' be a compact set. For Q <t < r. the (r, t)-doublc offset Kr^t of K 
is the set defined by: 

Kr^t = {KrYv 

Using a result of Federer ([H], Theorem 4.8) stating that the distance function 
to a closed set A with positive reach ri is differentiable with non-zero Lipschitz 
gradient on the complement of the closure of the medial axis of A, the following 
result is obtained in [7j: 



RR n° 6756 



8 



Chazal & Cohen-Steiner & Lieutier & Thibert 



Theorem 3 (Double offset theorem) If r < for some value /x > and if 
t < then dKr.t is a smooth C^'^-hypersurface. Moreover, 

reach(9Ki._t) > niin(t, fir — t) 

which implies that the smallest of the principal radii of curvature at any point 
of dKr,t is at least min{t, fir — t). 

3 Definitions and background on curvature mea- 
sures 

For more details on that section, one may refer to [TH [STJ [3T] . All the notions 
presented here are used in the proof of the main result (TheoremU]). Theorem^] 
gives a result of stability for the curvature measures. The reader only interested 
in the result stated in Theorem[4]can skip most of this section and go directly to 
the subscction l3.5l in that section, we give another way of defining the curvature 
measure that does not use the notion of currents and of normal cycles. 

3.1 General currents 

Let I?™ the R-vector space of C°° differential m-forms with compact support 
on M" (see for details on differential m-forms). can be endowed with 
a topology similar to the topology on the space of test functions used to define 
distributions, the so-called C°° topology. A sequence {(jii) of elements of 
converges to S I?™ in the C°° topology if and only if there exists a compact 
K containing the supports of all the (pi such that the derivatives of any order 
of the (t)i converge to the corresponding derivatives of ((> uniformly in K . 
The topological dual of is the R-vector space 2?m of m-currents on R". 
Equivalcntly, currents can be viewed as differential forms whose coefficients are 
distributions instead of smooth functions. The support Spt{T) of a current T 
can be defined as the union of the support of its coefficients. 

The subset of m-currents with compact support is denoted by f™. We endow 
'Dm with the weak topology; 

Urn Tj ^ T <^ ycj) e X>" hm ^,-(0) = 7(0) 

j — ^oo j — ^oo 

3.1.1 Operations on currents 

Basic notions relative to differential forms can be transposed to currents by 
dualization : 

1. Boundary : to each m-current T one can associate am— 1-current dT 
called its boundary, defined by : 

Vcf) e dT{(j>) = T{d(f>) 

where d(j) denotes the exterior derivative of </> (see [35] for the definition). 
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2. Push-forward : given an m-current T and a smooth map / defined on a 
neighborhood of the support of T, one can define the push-forward ffT of 
Thy f : 

Wcf, e fiTicf,) = T{f*cP) 

where f*4> is the puU-back oi (f> hy f (see [31]). Note that this definition 
only makes sense when f*(j) is compactly supported. We thus have to 
assume that / is proper, that is f~^{K) is compact for every compact 
K. Actually, for the push-forward to be defined, it is sufficient that the 
restriction of / to the support of T is proper. 

Since /* commutes with exterior differentiation, commutes with the boundary 
operator, so that the push-forward of a current without boundary is also without 
boundary. 



3.1.2 Current representable by integration 

We say that a current T € Vm is representable by integration if there is a Borel 
regular measure ||r|| on R" finite on compact subsets and a unit to- vector fields 

T defined almost everywhere such that 

V(/. e T((/)) = J < T,<?!) > d||r|| 

Currents representable by integration are analogous to distributions of order 
0. A current representable by integration T can be "restricted" to any ||T||- 
measurable set A (see pp 356). The obtained current T\_A is defined by : 

V0 e 23™ Tl^((/.) = / < r,(/)> 1^ d||T|| 



3.1.3 Rectifiable and integral currents 

In particular, one can associate an m-current representable by integration to 
any oriented TO-rectifiable subset S of dimension to of R" (see [H]). It is a 
well-known fact that rectifiable sets of dimension to have a well-defined tangent 
space at 7i™-almost every point. Let S be the unit to- vector field encoding 
these -oriented- tangent spaces. The current associated with S, still denoted by 
S, is defined by : 

= / <S,(I)> dH"" 
Js 

More general currents can be defined by incorporating integer multiplicities fi 
in the previous formula : 

T(0) = / ^ < ^, > dH" 
Js 

If the support of S is compact, and ^dTi.™ < oo, we say that T is rectifiable. 
The space of rectifiable currents is denoted by TZm ■ 

A current is said to be integral if it is rectifiable and if its boundary is 
rectifiable. 
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3.1.4 Mass and norms of currents 

The norm of a m-differential form (/) is the real number 

11011= sup ||0p||, 

where, for each p g M", 

||0p|| = sup{| < <j>pXp > Up e A™TpA./", \Cp\ = 1}. 

There are different interesting norms on the space of currents 2?^- We mention 
the main ones: 

• The mass of a current T S Vm is the real number 

M(r) = sup{r(0), such that e 2?"M|0|| < i.} 

For rectifiable currents of dimension m, the mass somehow generahzes 
the notion of m- volume : the mass of the current T defined by 13.1.31 is 
Jg ^dH"^. Rectifiable currents thus have finite mass. Using general results 
on representation theory of geometric measure theory, it can be proved 
that if M(T) < oo, T is representable by integration. 

• The flat norm of a current T G Vm is the real number 

T{T) = inf{M(A) + M{B) such that T = A + dB,Ae Um^B e Um+i}- 

It can be shown that the flat norm can also be expressed in the following 
way : 

T{T) = sup{r(0), such that e ,M \ < 1, |M</)|| < 1}. 
3.2 Normal cycle of geometric sets 

For many compact sets of R'', one can associate a (d-l)-current that general- 
izes the notion of unit normal bundle of a smooth manifold (see |24] for more 
details). In the case where C is a compact set of M'' enclosed by a smooth 
(d-l)-dimensional manifold, the normal cycle N{C) of C is just the integral (d- 
l)-current associated to its outer unit normal bundle: 

S{C) = {{p,n{p)), p e dC, n{p) is the outer unit normal at p}. 

If C is a convex set, the normal cycle N{C) is just the integral (d-l)-current 
associated to the oriented set: 

5(C) = {(p,n), pedC, neCNcip)}, 

where CNc{p) = {n € S^, Vx G C n.px < 0} is the normal cone of C at the 
point p. If C = U^i^iCi is a union of convex sets, its normal cycle is defined by 
additivity. By using the inclusion-exclusion formulae, we have: 

n 

7V(C) = ^(-l)'+i ^K=A). 

1 = 1 l<ii<...<i,<n 
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The last definition implies that the normal cycle is well defined for polyedron, 
but also for a finite union of balls. In fact, Joseph Fu [H] proved that the 
normal cycle can be generalized to a very broad class of objects, which he calls 
geometric sets. This class of objects contains in particular subanalytic sets |23) . 
definable sets [3], riemannian polyedra [57], or sets with positive reach 

3.3 Curvature measures 

Let us first recall some basic definitions and notations in the case where M is 
a smooth surface that is the boundary of a compact set V of M''. The unit 
normal vector at a point p € M pointing outward V will be refered as n{p). 
Note that M is thereby oriented. Given a vector v in the tangent space TpM 
to Af at p, the derivative of n in the direction u at p is orthogonal to n{p). 
The derivative DpU of n at p thus defines an endomorphism of TpM, known as 
the Weingarten endomorphism. The Weingarten endomorphism is symmetric. 
The associated quadratic form is called the second fundamental form. Eigen- 
vectors and eigenvalues of the Weingarten endomorphism are respectively called 
principal directions and principal curvatures. In the 3-dimensional case, both 
principal curvatures can be recovered from the trace and determinant of Dpn, 
also called mean and gaussian curvature. 

The curvature measures can be defined for any compact of R'' admitting a 
normal cycle (i.e. for geometric sets). Before introducing their definition, we 
first need to give the definition of invariant forms as follows. 

We identify the tangent bundle TW'' with E x F, where E is the base space 
and F is the fiber. Let J : E ^ F he the canonical isomorphism between E 
and F. We endow TR'^ with the dot product < (e, /), (e', /') >=< e,e' > + < 
J~^{f)J~^{f') >■ At any point (m,^) of STR'^ = {(m,e) e E x F, ||^|| = 
1}, we consider an orthonormal frame (ei, ...,ed-i) of the space orthogonal to 
J^^(^) and we take = J{ei). We build the (d-l)-differential form: 

r! = (e*+te*)A...A(e*_i+te*„J, 

where u* denotes the 1-form defined by u*{x) =< u,x >. One can show that 
this form does not depend on the chosen orthonormal frame. The coefficient of 
is a (d-l)-form denoted by uji. One can show that each tOi is invariant under 
the action of the orthogonal group. We now define the curvature measures as 
follows: 

Definition 3 Let C be a geometric compact subset of'R'^. The k*^^ -curvature 
measure ofC, denoted by <I>^ associates to each Borel subset B ofW^ the real 
number: 

=7V(C)l(BxF)K). 

Note that '^^{B) = N{C){1bxf i^^k), where Ibxf is the indicatrix function 
oi B X F. Therefore, we can extend this notion of curvature measure to any 
Lipschitz real function. This point of view is crucial and is a key point that will 
allow us to state simple results of stability in this paper. More precisely, one 
defines the isotropic curvature measure for every Lipschitz function / on by: 

$^(/) = 7V(C)(/c.,), 
where / is defined on R'* x R'' by f{p,n) = f{p). 
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One can show that if C is the volume enclosed by an hypcrsurfacc dC, then 
(B) is the integral over dCOB the k-th symmetric function of the principal 
curvatures of dC [31] . 

In dimension 3, and <I>^ are respectively the integral of twice the mean 
curvature and the integral of the gaussian curvature. We then take the notation 

= / G{p)dp and = / H{p)dp, 

J BndC JBndC 

where H[p) an Giji) denote respectively the gaussian and the mean curvature 
of dC at p. 

3.4 Anisotropic curvature measure 

The anisotropic curvature measure [HI dHl [31] is defined as follows: 

Definition 4 Let C he a geometric subset of . The anisotropic curvature 
measure ofC, denoted by He associates with every Borel subset B of the 
following bilinear form on E: 

VX,y e K'' Hc{B){X,Y) = N{Cy{B x F){uj^'^), 

where 

uj^'-"" = *e{J-\0* /\ X*) A j{Yy , 

here *e denotes the Hodge dual on E. 

In the particular case where C is the volume enclosed by a compact hyper- 
surface, Hc{B) is just the integral over i? n 9C of a symmetric bilinear form 
He related to the second fundamental form of dC. More precisely, this form He 
coincides with the second fundamental form of dC on the tangent plane of C 
and vanishes on its orthogonal complement. For any Borel set B of M'', one has 
[HIISIISI]: 

Hc{B) = f Hc{p)dp. 

J BndC 

In the 3-dimensional case, it is also convenient to introduce another form 
w"^'^ useful, in particular for polyedrons [15]. At a given point (m, ^) G K'^ x S'^, 
one has: 

^S) = (0' A (C X X, 0)* and ufj^^^ = (X, 0)* A (0, ^ x Y)\ 

where x is the cross-product in M?' . It is then possible to define another 

anisotropic curvature measure, denoted by Hq, that associates with every Borel 
subset B of M'^ the following bilinear form on E: 

VX,r gM'' Hc{B){X,Y)=N{Cy{B X F){oj^'^). 
In particular, if C is the volume enclosed by a compact surface, one has: 

Hc{B)^ I Hc{p)dp, 

J BndC 
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where He is defined as having the same eigenvectors than He , but with swapped 
eigenvalues on the tangent plane, and vanishes on the orthogonal component of 



Similarly as in the previous subsection, one can also define the anisotropic 
curvature measures Hc{f) and Hc{f) for every Lipschitz function / on WK 

3.5 Another approach for defining curvature measures 

It is possible to define the curvature measures of a set with positive reach as the 
limit of the curvature measures of its offsets. The advantage of these definitions 
is that they do not rely on the notion of normal cycle. 

More precisely, let y be a set with positive reach R> and let t < R. It is 
known that dVf is a,C^'^ hypersurface of R** [22]. The second fundamental form 
and the principal curvatures of dVt are thus defined almost everywhere. There 
is of course no pointwisc convergence of the principal curvatures when t tends 
to 0. However, the integrals of the curvatures of dVt converge to the integrals 
of the curvatures of V when t tends to (this is a consequence of Lemma H]) . 
That allows us to define the isotropic curvature measures of V for every Borel 
subset B of M."^ as follows: 



where pv is the projection onto V; B' = {p E W^, pv{p) £ B}; s'' is the 
k-th elementary symmetric polynomial of the principal curvatures Ai,...,Arf_i 
of dVt- In other words, s*^ satisfies for every a; G K: {x + Xi)...{x + Xd-i) = 
s° + s^x + ... + s'^^^x'^~^. Now, remark that we have: 



where Is is the indicatrix function of B. Therefore, we can extend this notion 
of curvature measure to any Lipschitz real function. This point of view is crucial 
and is a key point that will allow us to state simple results of stability in this 
paper. More precisely, one defines the isotropic curvature measure for every 
Lipschitz function / on by: 



Similarly, one extends the notion of anisotropic curvature measure of |14| ITS] : 
the anisotropic curvature measure of V associates to any Lipschitz function / 
the d X d symmetric matrix defined by: 



where Hgvt is a matrix-valued function defined on that coincides with the 
second fundamental form of dVt on the tangent space, and vanishes on the 
orthogonal component. 



C. 
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Now, let K be a compact set whose /i-reach is greater than r > 0. Then 
V = has a reach greater than pr. It is then possible to define the extended 
notions of curvature measures of by (see Lemma |4]) : 

= and HkM) = -HvU)- 

In the 3-dimensional case, it is also possible to define a second anisotropic mea- 
sure curvature H by: 

Hv{f)^lim f f{pv{p)) HqvAB) and 1 kM) = ^'HvU). 

where Hgvt is defined as having the same eigenvectors than Hgvt i but with 
swapped eigenvalues on the tangent plane, and vanishes on the orthogonal com- 
ponent of dVt. 

4 Stability results 

4.1 Curvature measures of the offsets 

The main contribution of this paper is Theorem HI Thanks to the formulation of 
the curvature measures with Lipschitz functions, the statement is simple. This 
theorem states that if two compact sets K and K' with positive /x-reaches are 
close in the Hausdorff sense, then the curvature measures of their offsets are 
close. We recall that the covering number M{A, t) of a compact set A is the 
minimal number of closed balls of radius t needed to cover A. 

Theorem 4 Let K and K' be two compact sets of whose ji-reaches are 
greater than r. We suppose that the Hausdorff distance e = duiK, K') between 

K and K' is less than min(/i, i). /// : R*^ M. is a Lipschitz function 

satisfying \ f\ < 1, then: 

\^'KM)-^K^if)\<kir,t,,dJ) sup(Lzp(/),l) Vi, 

and 

\\HKAf)-HK;Xf)\\ < kir,fi,dj) sup(Lip(/), 1) V^, 
where k{r, fi, d, f) only depends on f through the covering number Af {spt{f) q(^^-^ , fir/2); 
Lip{f) is the Lipschitz- constant of f ; spt{f) = {x S R"^, f{x) ^ 0}. 

The proof of this theorem is given in Sections El [7] and [HI We show in Figure 
[5] that this bound is tight. Furthermore, in the 3-dimensional case, this result 

also holds for the anisotropic curvature measure H. 

Now, if we take the function f{x) = max(l — — c||/r, 0) equal to 1 at 
a point c G dK'^ that radially decreases in a "small" ball B of radius r and 
vanishes out of B, then we can get local information about the curvature of K'^ 
from the curvature of Kr in the neighborhood of c. 

We also note that the conclusion of the theorem may be rephrased by saying 
that the bounded Lipschitz distance between the curvature measures of Kr and 
K'^ is bounded by 0{y/e). The bounded Lipschitz distance between measures 
is similar to the Wasserstein distance (also called earth's mover distance), ex- 
cept that it applies to general signed measures whereas Wasserstein distance is 
limited to probability measures. We refer to [S] for precise definitions. 
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Figure 2: Tightness of the bound: we take K = [p, q] and K' = [p, q]U{s}, where 
s is at a distance e from K. We have dH{K,K') = e and the total curvature 9 
of K'j. between a and b satisfies 9 — 2 arccos {^-^) = O {^/e) . 



4.2 General result 

The result of the previous section ensuring the stability of the curvature mea- 
sures assume that both the compact sets K and K' have sufficiently large /i- 
reach. Nevertheless, in practical settings, particularly when dealing with point 
clouds, such an hypothesis is rarely satisfied. Using Theorem [l] it is still possi- 
ble to approximate the curvature measures of the offsets of a compact set with 
positive /i-reach from any sufficiently close approximation of it. 

Theorem 5 Let K and K' be two compact subsets ofW^ such that r^{K') > r. 
Assume that the Hausdorff distance e = di{{K, K') between K and K' is such 

2 

that £ < go_^9^2 f- Then the conclusions of Theorem^ also hold. 

Of course, if for s > 0, the compact set K' satisfies rfj_{K'^) > c > 0, the same 
theorem applies to K'^. Furthermore, thanks to Theorem[l] the value of r^{K'^) 
can be read on the critical function of the sample K. 
Proof : 

It follows from Theorem[l](the critical values separation theorem of [6]) that the 
critical function of K is greater than ^ on the interval ^"^j r — S-^/ir^ . Note 
that since n > ^, the critical function of K is also greater than ^ on the same 
interval. As a consequence the two compact sets K — Kwe and K' ~ K'lee 

have their ^-reach greater than r — 'i\fer — Notice that for any S > 0, 
Kg ~ -f^'is'+i and since dniK, K') = e then dniK, K') < e. To apply Theorem 
m to and K' , ^ and r, the Hausdorff distance e between K and K' must 
satisfy < e < fi^ {r - 3^/er - ^) (note that since f < ^, min(/^, i) = /_*). 
One easily checks (by computing the solutions of the inequality and using that 
for all < .T < 1, VI ^ X > 1 — f ) that this inequality is implied by the as- 
sumption made on e. The theorem now follows immediately from Theorem!?] □ 

Remark that in the previous proof, the choice of ^ to apply the critical 
values separation theorem of [B] and then the theorem H] is arbitrary. The same 



RR n° 6756 



16 



Chazal & Cohen-Steiner & Lieutier & Thibert 



proof can be given replacing ^ by any < /i' < /i leading to somewhat different, 
but more technical, constraints on e and r in the statement of Theorem |H 

5 Computation of the curvature measures of 3D 
point clouds 

5.1 Gaussian and mean curvature measures of 3D point 
clouds 

Let K he a. finite set of points in and r > 0. We assume that the set of balls 
of radius r and centered in K are in general position (as in [17j): no 4 points 
lie on a common plane; no 5 points lie on a common sphere; and the smallest 
sphere through any 2, 3 or 4 points of K has a radius different from r. In the 
3-dimensional case, we denote by = the mean curvature measure and by 
~ tl^^ gaussian curvature measure. 

Note that the boundary of the union of balls Kr is a spherical polyhedron: 
its faces are spherical polygons contained in the spheres of radius r centered on 
points of K; its edges are arcs of circles contained in the intersection of pairs 
of spheres of radius r centered on K; its vertices belong to the intersection of 
three spheres of radius r and centered in K. It follows from Lemma 2.2 in [TH] 
that the combinatorial structure of dKr is in one-to-one correspondence with 
the boundary of the a-shape of K for a — r. As a consequence, under the 
general position assumption, to compute <i>^ (/) and (/), it is sufficient to 
compute the curvature measures ^^if) ^^'^ if)' where C is the union of one, 
two or three balls of radius r. The orientation of the boundary of the union of 
ball is taken so that the normal is pointing outside. All the curvature measures 
for the union of one, two or three balls are given respectively in Propositions [TJ 
[HandU 

Proposition 1 Let M be a ball of radius r o/M'^ and let B be a Borel set oJM?. 
Then the curvature measures of B above B are given by: 

2 1 

- Area{B n 9B) and ^^{B) = — Area{B n 9B). 

Proof : 

Since 9B is a smooth surface, one has 

^"{B) = / H{p)dp = - Area{B n 91) 

JBndB ^ 

and 

(S) = / G{p)dp = \ Area{B n dM), 

JBndB '"' 

where G{p) is the Gaussian curvature and H{p) is the mean curvature of B at 
p. □ 
In the following, we use the notation: ui^ = loi and w*^ = W2- 

Proposition 2 Let Mi and B2 be two intersecting balls ofM.^, of same radius 
r > and of centers Ai and A2, and C be the circle 9Bi n9B2. Let B be a ball 
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of M. . Then the curvature measures of Mi U B2 above B Ci C are given by: 



AiA2\ I fAiA2^^ 



KuM, n C) = -2/3r arcsin j y 1 y 2r 

and 

r , . P A1A2 



where (3 — 2ti /^p)^^ is the angle of the arc of circle B C\ 



C. 



Proof : 

Since the normal cycle is additive and since C is one-dimensional, one has: 

<uB, n c) = < (i? n c) + < (B n c) - $^^^1. n c) 
= -<nB.(snc). 

We now need to describe the normal cycle of Bi n B2 "above" C O B. Since 
Bi n B2 is convex, its normal cycle "above" C fl B is just the 2-current defined 
by integration over the set: 

^clBi nBa) = {(m,0, meCnB, \\£_\\ = 1 andVq e Bi n B2 tTiq.^ < 0}. 

Let a ~ arcsin (^^^^y In an suitable frame, the set 5'c(Bi n B2) can be 
parametrized by: 

/: [0,f3]x[-a,a] x 

cos u 



r cos a cos u \ 








r cos a sin u 


\ ,sin(u) j 


• 


+ cos(w) 












{u,v) 



Let (m, = /(w, v) e S'(Bi n B2). Wc put 

(— sin u \ ( ^ \ ( "^^^ ^' 

cosu and 62 = cos(?;) — sin(i;) sinu 
/ V 1 / V 

Then (ei,e2,^) is a direct orthonormal basis of M'^. We put t\ = (ei,0), £2 
(62, 0), tx = (0, ei) and = (0, 62). On has [l4] : 

= £1 A £2 + £1 A £2 and = £1 A £2. 
Furthermore, one has: 

df df 

— (m, u) = (r cos a 61, cos w 61) and — (w, u) = (0, 62). 

aw aw 

We then have: 

rc^<=((l,0),(0,l))=cosz;. 
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Then 

r0 pa 



^MiUMii^) — ^ / / r COS a dudv — —2(3r a COS a. 

Jo J -a 



Similarly, one has: 



*BiuB2 (C') = - / / COS u dudv = "2/3 sin a = 5—^. 

Jo J-a 

□ 



Proposition 3 Let Bi, B2 and B3 he three intersecting halls of , of radius 
r, of centers Ai, A2 and A^, and p £ dMi D dM2 H dM^. Then the curvature 
measures of Mi U B2 U B3 ahove p are given hy: 



and 



*BiuB2UB3({p}) = 4 arctan ^ tan tan — tan — tan 
where 



a — ai-j 



^ 01,2+02,3+03,1 



2J 



{pAi.pA^ = 2 arcsin (^^^|^) , 



Proof : 

Since the 2-form u)^ is mixed and since the support of the normal cycle of 
Bi U B2 U B3 "above" p lies in {p] x M^, one has $BiuB2({p}) = 0. Since the 
normal cycle is additive and since {p} is 0-dimensional, by using Proposition [51 
one has: 

*^iUB2UB3(M) = *^i(M) + *f2(M) + 'J'f3(M) 

-*|nB2(M) - *BinB3(W) - *f2nB3(M) 
+*^BinB2nB3({?'})- 
= '^BinB2nB3({^'})- 

*BinB2nB3(M) is just the area of the set: 

S{v) = {(p, 0, ll^ll = 1 o.nd Vq e Bi n B2 n B3 ^-C < 0} . 

The set Siji) is a spherical triangle whose area is given by the following formulae 
(see [I] page 289): 



4 arctan i/tan ( — | tan {'^ — \ tan {'^ ^^2^\ / Q^i.-i 



2/ \ 2 ) \ 2 I \ 2 



□ 



In Figure [31 we sampled a non-manifold compact set. It is interesting to 
note that both the mean and the gaussian curvature are low (and negative) in 
the middle of the disc. This is due to the fact that the curvature is "caught" by 
a function / whose support traverses the disc and intersects the cube. Figure [4] 



INRIA 



Stability of Curvature Measures 



19 




Figure 3: The mean (left) and gaussian (right) curvatures of a point cloud P sampling 
a non smooth compact set union of a solid cube with a disc and an arc circle. The 
colors vary from green (minimum value of the curvature) to red (maximum value of 
the curvature) passing through blue. The curvature measures are computed using our 
approach for a fixed offset value of q = 0.1 (the diameter of the point cloud being 
equal to 2). A curvature value is assigned to each data point p by integration of a 
Lipschitz function with support contained in a ball of center p and radius 0.3. These 
values are then used to color the faces of the boundary of the a-shape of P (top) or 
the offset of the point cloud (bottom). 

also illustrates the algorithm for a sampling of the union of a torus and a cube. 
The algorithm 

In practice, we use a halfedge data structure (provided by S. Loriot from INRIA 
Sophia- Antipolis) for the representation of the boundary of the union of balls Kr- 
Let now / be a Lipschitz function on Mf^. To compute the isotropic curvature 
measures $§-(/) and $^ (/) of the union of balls Kr, we triangulate dKr- 
More precisely, we triangulate every face of dKr into "small" spherical patches of 
diameter less than rj. We then approximate / by a function / that is constant on 
each "small" patch and that coincides with / in at least one point of each patch. 
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Figure 4: The mean (left) and gaussian (right) curvatures of a point cloud P sampling 
the union of a solid cube and a solid torus. The colors vary from green (minimum value 
of the curvature) to red (maximum value of the curvature) passing through blue. The 
curvature measures are computed using our approach for a fixed offset value of a = 0.1 
(the diameter of the point cloud being equal to 2). A curvature value is assigned to 
each data point p by integration of a Lipschitz function with support contained in a 
ball of center p and radius 0.3. These values are then used to color the faces of the 
boundary of the a-shape of P. 



Wc then calculate the isotropic curvature measures with / "above the faces" by 
using Proposition[T] Similarly, by using Proposition[2l we calculate the isotropic 
curvature measures above the "spherical edges" with a function / piecewise 
constant on "small" edges of lengths less than r]. Finally, we calculate exactly 
the isotropic curvature measures "above" the vertices by using Proposition [31 
Numerical integration 

Let B' be a borel subset of R"^ that contains the support of the functions / and 
let B = B' D dKr- The numerical error done by approximating / by a piecewise 
constant function / is then given by: 

\^kM) - 'fK/)! < Ltpif) M{N{Kry{B X R")) rj 

and 

- '^kM)\ < L^p{f) MiNiKrHB x K'')) r, 

where / is constant on each patch and coincides with / at at least one point in 
the patch. 

Proof: We show the result for the Gaussian curvature measure. The proof for 
the mean curvature measure is similar. One has: 

l*l(/)-<5'l(/)l = / 

J spt{N{Kr))n{BxR'i) 

< sup \\f{p,n)-f{p,n)\\M{N{KrHBxR'')). 

Let be a face of the spherical polygon dKr n B' , that is decomposed into 
spherical patches A. Let / be the piecewise constant function equal to /(po) on 
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every spherical patch A of F (with po E A). Wc then have: 

sup \\f{p,n) ~ f{p,n)\\ < Lipif) sup 77A, 

{p.n)e{FxR'') ACF 

where rjA is the length of the longest edge of A. The result is similar for a 
spherical edge C of dKr OB': if one decomposes C into small arc of circles e 
one has: 

sup \\f{p,n) - f{p,n)\\ < Lip{f) supZ(e), 

(p,n)e(CxR'') eCC 

where l{e) is the length of the arc of circle e. 

5.2 Anisotropic curvature measures of 3-D point clouds 

In the previous section, we described an algorithm that allows to calculate the 
curvature measures $^ (/) and <I>^ (/) when if is a finite set of points. Sim- 
ilarly, one can also calculate the anisotropic curvature measure HK^if) and 

Hx^if) of Kr by using the curvature measures formulae for the union of one, 
two or three balls (see Propositions [H [5] and E]) . More details on how to use 
these formulae can be found in [TH [H] . 

Proposition 4 Let M be a ball o/R^ and let B be a Borel set oflS?. Then the 
anisotropic curvature measures ofM above B are given by the 3 x 3-matrix: 

Hm{B) ^Hb{B)^ [ - {Id ~ Up' Hp) dp, 

JBndB ^ 

where Id is the identity matrix o/M"^, Up is the unitary normal o/B at p and 
^Up its transpose. 

Proof : 

Let X, Y be two vectors of R'^ and Pp denotes the projection onto the plane 
tangent to B at the point p. Using that Pp{X) = X— < X,np > Up, we have: 

Hb{B){X,Y)= [ Hp{Pp{X),Pp{Y))dp^ [ - {Id -Up' Up) dp. 

JBndB JBndB ^ 

□ 



Proposition 5 Let Mi and M2 be two intersecting balls o/R'^, of same radius 
r > and of centers Ai and A2, and C be the circle dMi n 9B2. Let B be a 
ball ofM.^. Then the anisotropic curvature measures 0/B1UB2 above BOC are 
given by the 3x3 matrix i?BiuB2(-S H C): 



I 



-r sma 



l-cos2ff 



1-cos 2/3 
4 



2 

sin 2q 











sin 2/j 



13 {a 



sin 2a 
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id the 3x3 matrix H^^yM^i^B n C): 



ra sm a 



I 



V 



sin 2g 
2 

cos 

2 





COS 2^-1 
2 

^_ sin2^ 







/ 



in the orthonormal frame (0, fc) where (0, generates the plane of the 
circle C and i points toward one on the two extremities of the arc of circle 
B n C, where a and {3 are defined as in the proof of Proposition\^ 

Proof : 

We take the same notations that the ones of Proposition [H Wc have: 

r0 



HM,uB,{BnC){X,Y) = 



f*LU 



X.Y 



As in Proposition H we take (m,^) = f{u,v) G S'c(Bi 082). Wc have: 



0), (0, 1)) = c J^^) (fi(«, v), ^{u, v) 



{0,Y) A X X, 0)((rsina ei,cos-(; ei), (0,62)) 
(0,y).(rsina ei,cosw ei) {£_ x X,0).(rsina ei,cos 

(0,y).(0,e2) (CxX,0).(0,e2) 
cosw Kei det(^,X, rsinaei) 
Y.e2 
= — rsina X.e2 Y.e2. 

The resuk is thus obtained by integrating the matrix of 62*62 given by: 



sin^ V cos^ u 



sm V sm u cos u 

— cos V sin V cos u 



sm V sm u cos u 



sin^ V sin^ u 



cos V sm V sm it 



' sm V cos w cos u 
- sinw sin u cos w 



Similarly, wc have: 



The result for i/j 
given by: 



((1,0), (0,1)) ==rsinQf X61 F.ei. 
JB2 (-Bn C) is thus obtained by integrating the matrix of 6i*ei 

sin^ u 
- sin u cos u 


□ 




Proposition 6 Let Bi, 182 and B3 he three intersecting balls of M.'^ , of radius 
r, of centers Ai, A2 and A3, and p e dRi D 9B2 dM^. Then the anisotropic 
curvature measures of Mi U B2 U B3 above p are equal to 0: 

^^BiUB2UB3 ({p}) = and iyBiUBgUBsClp}) = 

Proof : 

The two forms cj^'^ and ui-^'^ being mixed, they vanish on the normal compo- 
nent. □ 
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6 About the normal cycle of a compact of 

In this section, we give preliminary results for compact sets in M'': in Lemma 
m we bound the mass of the normal cycle of a set whose boundary is of class 
C^; in Lemma |3l we recall a result of Zahle et al. that allows to define, under 
several assumptions, the normal cycle of a set whose complement has positive 
reach; in Lemma UJ we give relationships between the normal cycle of a set 
whose complement has positive reach and the normal cycle of its offset. 

We first give the following lemma (that is a direct consequence of the Isotopy 
Lemma, Proposition 1.8 in [26]): 

Lemma 1 Let K be a compact of . If r is not a critical point of the distance 
function (i.e. Xk{'>') Oj, we then have: 



{.T, dxix) < r} = {x, dxix) < r} and {x, dxix) > r} ~ {x, dxix) > r}. 

Lemma 2 Let X be a compact subset of M'^ whose boundary is a hypersurface 
with positive reach R. We then have: 

d-1 



MiN{X))< [1 + — ] W'-'idX). 



More precisely, if B' is a Borel ofW^ and B ^ B' H dX, we have: 

d-1 

M[N{Xy{B xR'')] < (^1 + ' H'^-'{B). 

Proof : 

Let us consider the following map: 

f : B spt{N{X)) niB X R'*) 

771 i-^ (m,n(rn)) 

Since the reach, and therefore the curvature of the surface dX is bounded 
by R, it follows that the map m i— > n{m) is -^-Lipschitz. That implies that / 



is also Y 1 + -^-Lipschitz. Since the map / is onc-to-onc, the general coarea 
formulae (see Theorem 4.2 of [30] or Proposition 3.13 of [55]) states that: 



M{N{Xy{B X R'^)) = / dH''-' 

Jspt(N{X))n{BxS.''} 

= [ Jd-i{f) dn"-' 

Jb 



< / (1 + ]^) ' dn"-' 

d-1 

< (1 + -^') ' n"-'iB). 



B 



i?2 

where Jd-i{f) is the (d-l)-jacobian of the Lipschitz map / (see [5S], page 24- 
25). □ 
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Under certain assumptions, Rataj and Zahle [34j showed that the closure 
of the complement of a compact set with positive reach admits a normal cycle. 
More precisely, we denote by i the inversion map 

and by the push-forward by i. We define the erosions of X by X^^ = 
(^X'^ We then get the following lemma by combining results of Zahle et al. 
[331 El and [7]: 

Lemma 3 Let K he a compact set of W'' with fi-reach greater than r and let 
X — Kr ■ We can then define the normal cycle of X as the flat limit of its 
erosions as follows: 

T{N{X^,)-N{X)) -^,^oO. 
Furthermore, since X'^ has positive reach, it has a normal cycle and we have: 
N{X) = -i^N (X^ and M{N{X)) = M {N (X^)) . 

Proof : 

Proposition 2 and Theorem 3 of [34] imply that the normal cycle N{X) oi X is 
well defined and we have: 

T{NiX^,)~NiX)) -^,^0 0. 

Now, since X"^ has positive reach [7], Corollary 3.1 of (or Proposition 4 of 
[31]) imphes that: 

T [N (X^) - i-i^N ^,^0 0. 

We then have 

N{X) = -i^N (X^) . 

□ 

We now consider the following map 

Ft: M'' X R'^ W'^ X R"* 

{x,v) {x + tv,v) ' 

and denote by F^ the associated push-forward for currents. 

Lemma 4 Let V be a compact set of M."^ with positive reach R > 0. We then 
have for e < R: 

N{V) = F.,^ [NiV,)]. 

That result can be localized. Let B' be a Borel of and B = (B' x R'^) n 
spt{N{V)). We then have: 

N{V)^B = F.,^ [N{V,)^FZl{B)\, 
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In particular, if K is a compact set of with fi-reach greater than r, X = 
Kr, the set V = X'^ has a reach R > fir, and we have for e < R: 

N {X) = -z^N (X^ ^-H^-^i (^0] 
and M {N {X)) < {1 + e^)"^ M [N (X^ ,)] . 

Let now B' he a Borel ofW^ and B = (B' x W^) n spt{N{X)). We then have: 

N (X) lB = -tiF^^t [N (X^,) , 



d-1 



M{N{X)^B)<{l + e') ' M[N{X-,)^FI^{i{B))]. 
Proof : 

We notice that the restriction of F_c to the support of N {V^) is a one-to-one 
map: 

: spt{N{V,)) spt{N{V)) 
{x, v) {x — ev, v) 

The general coarea formula implies that: 

M{N{V))= [ dn''-^=( Jd-i{F_,) dw'-K 

JsntlN(V)) JsvtiNiV,)) 



d-l 



The map F_e is vT+e^-Lipschitz, thus we have Jd-i{F^f) < (l + e^) ^ and: 



M{N{V))<il + e^Y-'' M{N{V,)). 



spt{N{V,)) 
spt{N{V,)) 



More precisely, for any d — 1-differential form tp on x W^, one has: 
F_,,iV(ye)M =N{V,){F1,^) 

iFl,^Ue^,A...Ae^,_.,)dn'-\x) 

spt{N(VJ) 

^F^,i.){{DF.,{x)ef) A ... A {DF_,{x)e-,_,)) dW'-^x) 

spt{N{V,)) 

, w , f (jJf-.(x)ef)A...Apf-,(x)eg_i) 

J.-i(F..)(x) ^^_.(.) [\(^DF.,^x)ef) A ... MDF^^ix)e^,_,)\ 

J,„i(F_,)(x) ^^_,(,) (ef-^^^) A ... A ejlf ^) dU'-H^) 

(ef A ... A e^ -^) dH'^~'^{y) (by the general coarea formula) 

pt{N{V)) 

= NiV) (^), 

where A .. A e^_i is a unit (d-l)-vector associated with the oriented tangent 
space Tx {spt {N [V^))) and A ..A e^_-^ is a unit (d-l)-vector associated with 
the oriented tangent space Ty {spt {N {V))) and Jd-i{F-^)[x) = \{DF-f_{x)el) ^ 
... A {DF^^{x)e^_T^)\ 7i''~-'^-almost everywhere. That implies that: 

N{V)^F^,^N{V,). 

The local version is done the same way by noticing that the restriction of F-^ 
to the support of N {V^) \-FZ^{B) is a one-to-one map: 

F.,: spt{N{V,)^FI^iB)) ^ spi (iV (F) lB) 
{x, v) 1^ {x ~ ev, v) 

□ 



dH'^ 
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7 Normal cycles of the double offsets of two 
close compact sets 

The aim of this section is to prove the following proposition. It states that if 
two compact sets K and K' are close in the Hausdorff sense, then the normal 
cycles of their double offsets are close for the flat norm. This result is local since 
we evaluate the normal cycles locally. 

More precisely, let B' be a Borel subset of and B = B' n dK'^ ^. We 
consider in the following the normal cycles of Kr^t and K'^ ^ restricted to B 
and tt{B), where tt is the projection onto ifr.t- We consider the two currents 
D' = A^(/C;_j)l(B X M'^) and D = N{Kr,ty{7r{B) x M'') given by: 

D'{^) = [ ^ 
and D{(p) = ip, 

■/spt(7V(ifr,t))n(7r(B)xE<i) 

for any (d-l)-diffcrential form of K'' x M''. 

Proposition 7 Let K and K' be two compact sets of M."^ whose ji-reaches are 
greater than r. We suppose that the Hausdorff distance e = dniK, K') between 

K and K' is less than ''^ min(/i, For t G 

can write: 

D' - D = dRi + i?2 
where Ri and i?2 are currents satisfying 

M{Ri) < M{D')kik2 
A/(i?2) < M{dD')kik2, 

where 

In particular the flat norm of D — D' is bounded by (M(D') + M(dD'))kik2. 
The proof of that proposition relies on the two following lemmas: 

Lemma 5 Let X and X' be two compact sets of R'^ with reaches greater than 
R > 0. Let t < ^ and e = dH{X,X'). If e < |, then for any x at a distance t 
from X , we have: 

2sm ^ <30^-. 



rji 



min(Aii) (2-\/2)' 2 
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Proof : 

For p such that < p < t, wc denote by Gx{x,p), as in section 5 of [5], the 
convex hull for every y g B{x,p) of all the "classical" gradients Vx(y) of the 
distance function dx- By using Theorem 5.6 in [5], we have that: 

Vx' (x) eGxix, p) p ■ (3) 



On another hand we know ([22] page 435) that the projection map ttx on X 
R 

R-(t+p)' 

\ rwid Vioc rlfnT- 1 1\ n on/^ rl f 'TT ^ f nf\ Tr-,^(niW : 

(*+p) 

Rp 



is j^_(^^p) -Lipschitz for points at distance less than {t + p) from X . Then, for 
y e B{x,p) one has d[x,y) < p and d{Trx{x),'!Tx{y)) < R^lt+p) 



Wv^xiv) - xt:x{x)\\ < p 



R-{t + p)' 



Using the fact that Vx(^) = "'" > zTTxjz) one get, for y S B{x,p): 
||Vx(2/)-Vx(x)|| < ^ ' 



t - p V R-it + p) 

This and Equation® gives: 

2sin ±^±^L±± = \\Vx{x)-Vx'{x)\\< ^ 



Taking p = \JH and e < | one gets: 
2sin^X^M^^< 




□ 



Lemma 6 Let K and K' he two compact sets ofW^ whose p-reaches are greater 
than r. Then for every t £ (0, rp), we have: 

dn {Kr,t , K t) < - and dn {dKr,t , dKl ^) < - . 

p p 

Proof : 

First remark that if we take two compact sets A and B with /i-reach greater 
than r, such that dniA, B) < e, one has: 

dH {Wr,iB;r) < ^ (4) 

Indeed, let x € Ar'^ . Then d{x,A) > r and d{x,B) > r — e. Let s i-^ o'(s), 
(t(0) = a; be the trajectory of Vb issued from x and parametrized by arc-length. 
While a{s) G Br, we have 

dB{a{s)) = deix) + f \\V B{(T{s))\\ds >r-i+sp. 
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Wc then have a{s) e {BrY for s > 4. As a consequence, there exists x' G {B^Y 
such that d{x,x') < i. We apply Equation ^ with A = K and B ~ K' and 
we get: 

dH{p<:Y.jKly)<-^ (5) 

We apply again Equation (j4|) with A = [K^Y and B ~ {K'rY '^ith e = - and 

Remark also that for any compact sets A and B, one has dniAt, Bt) < dniA, B). 
Therefore, by Equation ([5]) one has: 

dH{Kr^t,Kt) < -■ 

The two last equations imply that 

dHidKr,t,dK^t) < -■ 



Indeed, let x € dK^.t- Then a; is at a distance less than ^ from K'^^ and K'^Y- 
Then, there exists y g K'^ ^^ and z £ K'^ Y such that xy < and xz < j^. Since 
the line-segment [yz] intersects dKr,t, there exists x' € dK'^ ^ such that xx' < -. 

Proof of Proposition [7] 

Let K and i^' be two compact sets whose /x-reaches are greater than r. Wc put 
e = dniK, K'). We suppose that — ^ ^ < t < Lemma [5] then implies 

that 

dH{dKr,u dKl t) < - < (2 - ^/2) 

Let f7t be the tubular neighborhood of dKr.t of radius t < reach{dKr^t)- The 
projection map tt onto ^-ftTr.t is then defined on Ut- We clearly have dK'^ ^ C Ut- 
More precisely, the map tt induces a one-to-one map between dK'^ ^ and dKr^t 
(see Theorem 4.1 of [H]). We now define 

tA: C/txM'' ^ spt{N{Kr^t)) 
(x,n) 1-^ (7r(a:),n^(^)) 

Let h be the affine homotopy between ip and the identity (see Federer [21] 4.1.9 
page 364) 

h : [0, 1] X {Ut X R'^) ^ R'^ X M'' 

{t,x) ^ (1 - t).T + t^/>(a;) ■ 

Since the map tt : ^ — > dK^^t is one-to-one, the map "0 also induces 
a one-to-one map between spt{D') and spt{D). Therefore, similarly as in the 
proof of Lemma m one has: 

4>^D' = D. 

According to Federer ([12], 4.1.9 page 363-364), we have: 
D-D' = YiD' - id^D' = dRi + R2 
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where Ri = [/itt ([0, 1] X D')] and i?2 hfi[0, 1] x dD'). Agam by Federer ([22], 
4.1.9 page 364) or Fanghua ([501 page 187), we get: 

Mih^{[0,l]x D')) <AIiD') sup sup |l + Jd„i(^^)|, 

spt{D') spt{D') 

and 

M {m[0,l] X dD')) <M{dD') sup jV" - sup |1 + Jd_i(V')|, 

spt{D') spt{D') 

where Jd-i{ip) stands for the (d-l)-dimensional jacobian (as defined in [29] . 
page 24-25). By Lemma HI the space component oi ip — id is less than ^. By 

Lemma [SJ the normal component of ip — id is less than 30. Thus 



sup \ip — id\ < 

spt{N{K'^^)) 



900- 



fit 



We can note that the jacobian of the space component of tJj (i.e. the projection) 



is bounded 



The jacobian of the map x e Kr^t is upper 



bounded by (7)'' ^ (where i < ^ is the reach of dKr,t)- The jacobian of the 

d-1 

We 



then have 



normal component of ip is upper bounded by the product t^^) 

\i + jd-im < 1 + 



i+t^ 



sup 

spt{N{K',^_,)) 



By using the definition of the flat norm we have: 

T{D-D') < Mih^ ([0,1] X D')) + M {h^{[Q,l] x dD')) . 
which is upper bounded by: 



/ 



{M {D') + M {dD')) 



I 



l + <^ 



t-^ 

A' 



V V 

8 From double offset to offset 



900- 



lit 



The aim of that section is to prove Theorem [H We assume that / is , the 
general case being obtained by a limit argument. Let C = {spt{f)a x R'^) n 
spt{N{K[)), where a is a positive real number that will be specified later and 
let B ~ Ft o i{C). Thanks to LemmajH we have: 

N {!<',) ^{zoF^tiB)) ^-tiF.ti[N{Kt)^B] . 

and 

N{Kr)^{ioF^tWB))) =-iiF^ti[N{Kr^t)^4'{B)]. , 
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where V' is defined in the proof of Proposition!?] Remark that we have N j) i_B — 
N {K'^j) l(B X R''), and N {Kr^)^^{B) = N (Kr.t) ^{^(13) x M^) where B is 
the projection of B onto the first factor of R'^ x and tt is the projection onto 
dKr^t- Hence, Proposition [7] aUows us to write: 

N (Kl) l(i o F^tiB)) - N (Kr) l(i o F^ti^iB))) = -i^F^^{dRi + R2), 

where Ri and R2 satisfy certain properties. In particular, from the proof we 
see that R2 = h'^[{N{K[. f)^dB) x [0,1]], where h is the result of the linear 
interpolation between the identity and -0. Since i o F-t is linear, we have 

[i o F^t) o /i = /i' o (i o i^„t, /rf[o,i]) , 

where h' is the result of the linear interpolation between the identity and ?/;' — 
ioF^toipoFtoi. We then have that -^F^^{R2) = h'^[{-i^F^tiN{Kl.j.),_dB) x 
[0,1]]. Noting that -^F_tiiN{Kl,tydB) = N{K'^ydC, and choosing a > 
\W — Id.\\oD, we see that the support of i^F-tiR2 is included in the a-offset of 
spt{N {K'^\dC) . Due to the choice of a, this subset does not meet spt{f) x R"*. 
As a consequence, for any differential form cj, letting: 

A = \<N (Kl) l(* o F^t{B))Ju >-<N (Kr) l(* o F^t{ij{B))), fuj > |, 

we have: 

A < \ < iiF^ttdRiJiu > \ 

< \ < iiF_tiRi,df ALU + fduj > \ 

< M{tiF_tiR,){Lip{f)\\u^\\ + \\du;\\) 

< {l+tY/'M{R,){Up{f)\M + \\du;\\) 

< M{N{K^,yB)hk2il + tY/\Lzp{f)\\Lo\\ + \\du\\), 

the last two inequalities following from Lemma HI Proposition [3 and the fact 
that Lip{f ) = Lip{f). Now remark that i o F^t{B) = C and i o F^t{i'{B)) = 
i]'{C) respectively contain spt{N{K';))C^(spt{f)y.W^) and spt{N{Kr))n{spt{f)x 
M''), due to the choice of a. Hence we have: 

A = \<N (/<,) ,fu>-<N (Kr) , fio > |. 

The quantity M{N{K'^ f)LB) may be bounded above using Lemma[2l 

M{N{KltyB) < {l + reach{dKl.ty^)'^n'^-\{dK',^t)<^B'), 

where B' is the preimage of spt{f)a H K'^. under the projection on K'^. Since the 
reach of dK^.^ may be bounded by a function of r and /i for t ~ rii/2, it only 
remains to bound the area term to finish the proof of Theorem U) To do so, we 
use (a local version of) Theorem H.3 in [S], which yields the following bound in 
terms of covering numbers: 

H''-\{dK^t) n B') < N{spt{f)^ n dK,t) ujd-i{2t) 

where ujd-i{2t) denotes the surface area of the ball of radius 2t in K''. Note 
now that M{spt{f)a n dK'^, t) < J\f{spt{f)a,t). The fact that a = 0{^/e) easily 
follows from the following lemma: 
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Lemma 7 

-Id\\^ = 0{V€). 

Proof : 

Let {y', n') e C and x' G dK'^ ^ be such that the projection of x' onto K''^ is y' . 
We put X — tt{x') the projection of x' onto dKr^t and we denote by n the normal 
to Kr^t at the point x that is pointing inside. We denote by y the projection of 
X onto K^. We then have: 

{x' , ~n') ~ Ftoi[y\ n'), [x, —n) = 4'{x' i ^n') and (y, n) = ioF^t{x, n) = ip'{y\ n'). 

We now need to bound yy' and ||7i— 7i'||. Lemnia[5]imphes that ||n— 7i'|| = Q{^/e). 
We introduce y the projection of x' onto if^. Let us now bound y'y and yy. 



y' 

Figure 5: Proof of Lemma [7] 
Since x' G dK'^.^, we have x'y' = i. Furthermore, by using Lemma [6l one 



has 5 — x'y G t — ^, i + ^ . Let now denote by y" the point of the half line 

issued from x' and passing through y', that is at a distance 5 from x'. We have 
by Lemma m 

y'y < y'y" + y"y 



< 



- + 2 sm — ^ ^ {t + - ] 

< - + 30 , A- (t+-) 

Now by using again Lemma El one has xx' < — and also that x and x' are at 
a distance less than (t+ - ) from K^. Since ilT^ has a reach greater than ^r, 



then the projection map onto is - — ^^^j-^j-Lipschitz for points at a distance 
less than (t + from (g2j page 435). We then have: 

e fir 

yy< 7 

M f.r-{t + f 

We then clearly have yy' < y'y + yy = O(v'e). □ 
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10 Conclusion and future works 

We have introduced the first notion of anisotropic curvature measure which is 
Hausdorff stable and applies to a large class of objects, including non manifold 
and non smooth sets as well as point clouds. Indeed, it is enough to require that 
some offset has a positive /i-reach, or, equivalently that the critical function of 
the set is greater than some positive number ^ on some interval. 

In light of these results, one can introduce a scale dependent variant of the 
normal cycle. We say that a compact set K C'R'^ satisfies the (Pa) property if its 
critical function is greater than some positive fj, on an open interval containing 
a > 0. For such K, we define the a-normal cycle as the rectifiable {d — 1)- 
current of x E'^-'^, N^iK) = F_atN{Ka). The effect of push-forward F_„tt 
is to move the support of the normal cycle closer to K: in simple cases (but 
not in general) Na{K) is equal to the normal cycle of the double offset i^a.a"^- 
This current captures in some sense the curvature information at scale a and 
has two nice properties. First, it coincides with the usual normal cycle for sets 
with positive reach, more precisely, if a compact set K has a reach greater than 
a, then Na{K) = N{K). Second it is Hausdorff stable, more precisely, if K 
satisfies {Pa), then there are constants C and eo > depending only on K 
such that if K' is a compact set such that dH{K,K') < e < cq then Na{K) 
and Na{K') differ by less than Cy/e in the so-called flat norm (see for example 
[551 HO] for a definition) , which implies that the associated curvature measures 
are also O {y/e) close. 

We think of several possible future directions. First we think of extend- 
ing our paradigm to the measure of higher order quantities such as torsion of 
curves or curvature derivatives. A possible track to define a stable measure for 
these higher order quantities is to integrate the gradient of a Gaussian function 
again the normal cycle: this would retrieve informations about the gradient of 
curvature measure. Associated stability results require more investigations. 

Furthermore, although the practical setting is beyond the scope of this paper, 
we obtain promising results for the estimation of the curvature measures from 
a noisy point cloud sample and we expect potential applications in the context 
of point cloud modeling. 
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